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ABSTRACT 

•  (Distribution  Limitation  Statement  No.  2) 

Charge  axchanga  cross  sections  for  fast  heavy  Ions  representative  of 

*  several  groups  In  the  periodic  table  with  atmospheric  and  other  atoms 
have  been  calculated  by  using  a  quanta  1  two-state  one-electron  method.  The 
principal  ions  considered  have  been  singly  charged  K,  I.  and  Cr  while 
target  atoms  and/or  amleculea  have  been  N,  0,  Ar,  Ne,  H  and  He.  Singly 
charged  lithium,  sodium,  rubidium,  and  cesium  have  also  been  considered 

as  projectiles  in  the  target  gases. 

Stripping  cross  sections  for  fast  heavy  Ions  incident  on  atmospheric 
and  other  atoms  have  been  calculated  by  using  an  Independent  scattering 
awdel  based  on  the  Born  approximation . 

The  relationship  between  the  quantal  predictions  snd  available  exper¬ 
imental  results  la  Illustrated  and  certain  improvements  to  the  simple 
theory  are  discussed.  A  comparison  la  also  made  between  the  electron  loss 
theory  and  available  experimental  results. 

Capture  Into  excited  states  la  discussed, and  estimates  are  made  of  the 
appropriate  cross  sections. 

9 

The  rate  of  energy  loss  by  heavy  ions  In  air  is  discussed. 
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SECTION  I 


INTRODUCTION 


In  this  work,  a  model  explored  In  reference  1  he*  been  applied  to 
the  calculation  of  further  charge  exchange  croaa  aectiona  for  colllalona 
between  heavy  Iona  and  atmoapheric  Iona  and  moleculea.  Modifications  to 
the  aimple  theory  are  auggeated  and  dlacuaaed.  The  method  uaed  la  the 
la  model  two-state  quantal  description  applied  to  the  ground  atatea. 
Estimates  are  alao  made  of  the  capture  cross  sections  into  excited  states 
of  the  projectile  atoms. 

Calculations  of  the  cross  section  for  electron  loss  have  been  per¬ 
formed  by  using  a  modification  of  a  theory  developed  by  Nikolaev  and 
Dmitriev  (Ref.  2).  In  this  model  the  projectile  electrons  are  assumed 
to  be  each  scattered  independently  by  the  target  atom,  according  to  the 
Born  approximation.  In  reference  2,  a  final  integration  over  a  simple 
analytic  form  factor,  which  is  replaced  herein  by  numerical  quadrature 
over  any  tabulated  form  factor.  Is  required. 

The  theory  of  charge  exchange  and  refinements  thereto  are  discussed 
In  Section  II.  Section  III  contains  a  discussion  of  the  electron  loss 
theory.  In  Section  IV,  available  experimental  results  are  compared  with 
the  calculated  results  while  the  Appendixes  describe  the  computer  programs 
developed . 
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SECTION  II 


QUANTA!.  THEORY  OF  CHARGE  EXCHANGE 


1.  Theory  of  Fast  Collisions 

A  brief  review  !•  preeenCed  of  Che  theory  given  In  reference  3 
as  Applied  Co  a  Cwo-sCeCe  one-elecCron  system.  Atomic  unlCs  are  em¬ 
ployed  ChroughouC.  Consider  point  A  Co  move  pasC  polnc  B  ac  a  conscanc 
velocity  v,l.e.,  a  rectilinear  trajectory  characterized  by  an  impact  param¬ 
eter  p.  Lee  a  single  electron  have  position  vectors  ra,  £b  and  £  with 
respect  to  A,  B  and  the  midpoint  of  AB,  respectively.  Let  the  electron 
see  potentials  Va(ra)  and  V^fr],)  centered  on  the  appropriate  polnts,and 
let  there  exist  bound  state  wavefunctions 

[t  +  Vra)  -  1  ^(r.)  «XP  <-  l  Et  t)  -  0  (1) 

considering  A  to  be  at  rest  and 

[T  +  vb<rb)  -  1  ^1  (Pj<£b)  **p  <-  1  Ej  t)  -  0  'Z, 

considering  B  to  be  at  rest  where  T  is  the  kinetic  energy  operator.  If  the 
midpoint  of  AB  Is  taken  to  be  at  rest,  wavefunctions  simply  derived  from 
<pt  and  cpj  are  solutions  of  the  time  dependent  Schrodinger  equation;  thus: 

[t  +  Va(ra)  -  l  r>  exp[-l(El  +  i  v2 )  t  ]  -  0  (3) 

and 

|T  +  W  -  i  £]  V£b-  expl'i(EJ  +  I  v2)cI  '  0  (4) 

where 

<?!<£».  O  “  'Pi<£«)  exp("  J  1  2  •  x)  <5> 

and 

cpjCEb,  JE>  “  'VjCEb)  e*P  (j  1  X  '  *)  <6) 
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The  total  wave  function  T  in  this  Approximation  Is  then  expanded  as 


*  -  a1<Pl  exp | -1  +  y  v2  )  t]  +  exp[-l  (ej  +  £  v2  )  tj  (7) 

and  substituted  Into  the  conditions 

[v  (« - 1  i)  »l«p[l(Ei  +  ?  v2)t|  -  o  (8) 

and 

h  *  (H  - 1  £)*1  “»Ii(ej  ■  f  v2)cl  -  °  <9> 

where  H  ■  T  +  V,  +  Vj,  Is  the  total  Hamiltonian.  The  equations  of  the  two- 
state  model  then  follow 

1  *i  +  1  StJ  exp[l(E1~  “  \t*t  +  htJ  exp[l(E^-  E^tjbj 

1  SJt  exp[l<EJ-E1)t]il  +  1  bj  -  hJt  expfKEj-E^tJaj^  +  h^bj  (10) 

where 

StJ  -  (<Pt.  <Pj).  htJ  -  (<plt  Va  (Pj),  and  h^  -  Cp^  Vb  cpt)  (11) 

The  Initial  conditions  for  Equation  (1)  eve  that,  say,  at  -  1,  bj  »  0.  At 
the  eod  of  the  trajectory  |bj|^  will  be  tiie  capture  probability  for  that 
Impact  paraswt'r. 

c(p)  -  |bJ(+«)|2  (12) 

and  the  cross  section  Is 


Q  -  2 n  m  c(p)p  dp  (13) 

Jo 


Normalisation  (-  |at|2  +  |bj|2  )  will  remain  at  unity  If 

hji  *  hu  “  1  *  <Ei  -  Vsi j  <14> 

which  replaces  Equation  (12)  of  reference  1. 
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Equation*  (10)  are  equivalent  to  the  following  form: 


-  1  Klj  b' 


exp^^*  8(t')dt'jexp(l(E1  -  Ej)t] 
bJ  “  "  1  Kjl  *1  J*  Stt'jdt'JexpCKEj  -  Et)t] 


where 


„  hk/  ~  Sk/  h£t 

^  •  -  lv,l* 


(15) 


(16) 


and  the  distortion  term  1* 


hil  ~  hJJ  ~  S1J  hJl  +  SH  hil 


1  -  S 


(17) 


i)1 


The  wavefunctlon  model*  and  the  appropriate  techniques  for  calculating 
the  matrix  element*  (11)  have  been  reported  in  detail  in  reference  1.  The 
method  of  integrating  Equations  (10),  the  set  solved  in  practice,  is  also 
described  in  reference  1.  The  wavefunctlon  parameter  in  the  Is  model 
employed  is  the  present  set  of  calculations,  and  others  available  for 
further  calculations  are  listed  in  Table  I. 

TABLE  I 


IONIZATION  POTENTIALS  AND  Is  MODEL  PARAMETERS 


Element 

Ionization  Potential  (au) 

a 

2a3/2 

Ki 

0.22009 

0.66347 

1.08083 

Fe 

0.29039 

0.76209 

1.33058 

U 

0.22463 

0.67027 

1.09750 

N 

0.53435 

1.03378 

2 . 102 19 

0 

0.50072 

1.00072 

2 . 002 15 

Ar 

0.57946 

1.07653 

2 .23390 
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2 .  Exchange  and  Distortion  Refinements 


Two  methods  are  described  below  by  which  it  is  believed  that  the  basic 
Is  or  ns  model  could  be  improved  should  further  calculations  of  this  type 
be  undertaken. 

a.  Electron  Exchange 

Let  u  represent  certain  matrix  elements  that  occur  in  Equations 
(10)  (u  m  Sj4,  hij,  or  h<i).  For  convenience,  the  complex  exponentials  in 
Equations  (10)  are  Included  in  the  definition  of  the  elements  in  this  anal¬ 
ysis.  The  various  approximations  to  p  are  denoted  as  follows: 

M  :  matrix  elements  calculated  by  using  exact  many-electron  fully 
entisymmetr iced  combinations  of  atomic  wavefunctions . 

M°  •  matrix  elements  calculated  by  neglecting  atom-atom  electron  ex¬ 
change  but  otherwise  exact. 

p  :  matrix  elements  adopted  for  Inclusion  in  the  two-state 
calculation. 

:  matrix  elements  calculated  directly  from  the  one-electron 
model . 

It  was  established  in  an  earlier  final  report  (Ref.  4)  that  at  suffi¬ 
ciently  large  lnternuclear  distances  where  second  and  higher  order  overlap 
terms  may  be  neglected 

M  -  F  u°  (18) 


where  F  is  the  geometric  mean  of  the  numbers  of  active  electrons  in  the 
initial  state  of  the  target  and  the  final  state  of  the  projectile.  An 
active  electron  could  be  defined  arbitrarily  as  one  with  a  binding  energy 
not  more  than  twice  that  of  the  most  loosely  bound  electron.  Because  the 
matrix  elements  u°  are  an  approximation  to  the  u°,we  might  write  analo¬ 
gously  to  Equation  (18),  as  an  approximation: 

U  “  F  u°  (19) 

but  at  the  smaller  lnternuclear  distance  F  u°  will  not  represent  u  well 
because  of  the  overlap  approximation  Inherent  in  F  and  because  of  the 
shortcomings  of  the  one-electron  model  Itself.  Therefore,  Equation  (19) 
is  modified  by  inserting  an  additional  factor  G  and  adding  a  correction,  & 
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(20) 


ii  -  G  F  40  +  6 


where  G  departs  from  unity  only  to  prevent  (I  from  becoming  unphyslcal 

(i.t.,  when  |F  j |  exceeds  unity).  The  correction  &  ensures  that 
normalization  is  conserved  so  that  [compare  Equation  (14)]: 


The  u°  elements  automatically  satisfy: 


d  S 


11 


dt 


(21) 


(22) 


and  Equation  (20)  is  applied  by  assuming 


SU  "  G  F  *1J 


hlJ-GFhiJ+2  A 


(2  3) 


VgfI5-?a* 


Substitution  of  Equations  (23)and  (22)  into  Equation  (21)  yields: 

A  “  7  if  F  (24) 

Now  the  limit  1b  imposed 

|s7jl  (25) 

where  L  might  be,  say,  0.9.  Then,  either 
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G 


1;  A  -  0 


(26) 


or 


G 


A 


i 


L  S° 

|S°| 


(27) 


This  extension  of  the  basic  model  contains  the  arbitrary  limit  L  and  the 
validity  of  the  scheme  requires  demonstrated  Insensitivity  of  the  cross 
section  to  L.  clearly, at  sufficiently  large  velocities  so  that  Equation 
(26)  holds,  the  effect  of  the  scheme  is  simply  to  multiply  the  cross 
section  by  F2. 

b.  Distortion 

At  large  internuclear  distances , the  rate  at  which  the  distortion  g 
diminishes  is  controlled  by  hj_j_  -  h j j  [see  Equation  (17)].  The  basic  one- 
electron  model  has  the  flaw  that  h^^  -  hjj  ~  (Zg-  Z^)/R,  assuming  singly 
charged  ions  in  both  channels,  whereas  in  fact  h^  -  hjj  should  diminish 
faster  than  R-*.  In  the  hope  of  improving  the  accuracy  of  the  results  at 
low  impact  energies,  the  program  has  the  capability  of  calculating  h^]_  and 
hji  with  the  Coulomblc  term  of  each  potential  set  equal  to  the  net  change 
of  the  appropriate  positive  ion.  In  the  calculation  of  h^j  and  hji  the 
original  form  of  the  potential  is  retained.  This  has  been  applied  to  + 
Ar  and  has  the  effect  of  reducing  the  cross  sections  somewhat  at  all  energies 
however,  this  does  not  constitute  a  clear  improvement  in  the  results  of  the 
calculation.  It  has  been  established  in  other  papers  (Ref.  5  and  6)  that 
distortion  has  a  great  effect  on  the  cross  section  at  lower  energies. 


3 .  Capture  into  Excited  States 

The  rate  of  energy  loss  of  a  beam  of  heavy  particles  in  gases  is  de¬ 
pendent  on  the  cross  sections  for  the  numerous  inelastic  scattering 
processes  that  can  occur.  The  contribution  to  the  rate  of  energy  loss 
from  charge  transfer  processes  into  excited  states  is  different  than  that 
for  the  capture  into  the  ground  state  of  the  projectile  system.  Capture 
into  excited  states  can  produce  excited  or  metastable  species  that  be¬ 
have  sufficiently  different  than  ground-state  systems  in  subsequent  pro¬ 
cesses  so  as  to  complicate  the  interpretation  of  many  experimental  studies 
of  heavy-body  collision  processes.  In  addition,  though  capture  into  excited 
states  and  possible  excitations  of  the  target  system,  energy  can  appear  in 
the  radiation  field  rather  than  in  the  potential  and  kinetic  energy  of 
the  nuclei  and  electrons.  The  majority  of  experimental  studies  of  charge 
transfer  processes  give  information  on  the  sum  of  the  cross  sections  for 
transfer  into  all  of  the  bound  projectile  states  with  the  target  left  in 
the  ground  state  or  excited  states,  while  the  theoretical  studies  usually 
consider  specific  states  of  the  final  target  and  projectile  system. 
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A  considerable  body  of  theoretical  and  experimental  literature  on 
charge  transfer  cross  sections  for  capture  Into  excited  states  is  avail¬ 
able,  however  it  Is  almost  exclusively  concerned  with  proton  projectiles 
on  various  target  gases,  so  that  generalizations  to  more  complex  projec¬ 
tile  systems  must  be  performed  with  care.  The  large  amount  of  data  for 
the  case  of  the  proton  projectile  indicates  quite  clearly,  however,  that 
the  cross  section  for  capture  into  an  excited  state  Is  usually  much 
smaller  than  the  cross  section  for  capture  in  the  ground  state,  and  that 
the  dominant  contribution  to  the  total  capture  cross  section  usually 
arises  from  capture  into  the  ground  state. 

If  the  cross  section  for  capture  from  a  state  with  principle  quantum 
number  n^  and  orbital  angular  momentum  quantum  number  Ij.  into  a  state  with 
quantum  numbers  nf  and  ff  is  denoted  as  Q  (n^  |n.if  ) ,  Oppenhelmer  (Ref.  7) 
showed  that  Q(ls|n£s)  falls  off  as  nf-3  at  hlghenergles  by  using  the 
Oppenhe imer -Brinkman -Krame r s  (Ref.  7  and  8)approxlmation, which  neglects  the 
nucleon-nucleon  interaction  in  a  Born-type  appoximation.  Generalizing 
this.  Bates  and  Dalgarno  (Ref.  9)  and  Bates  and  McCarroll  (Ref.  3)  argue 
that  Q(n£s|nfs)  falls  off  as  (nfn^)-3  at  high  energies  and  Q(n£l£|n£l£) 
for  /i,/f  4  0  falls  off  rather  more  slowly.  Mapleton  (Ref.  10),  in  full 
Born  approximation  calculations  for  protons  on  atomic  hydrogen  and  helium 
obtains  numerical  results  that  confirm  the  rules  given  above.  In  an 
experiment  Involving 

H+  +  He  -  H(3s)  +  He+ 


Hughes,  et  al.(Ref.  11)  obtain  results  that  approach  Mapleton's  results 
at  high  energy,  but  are  much  lower  at  low  energies  as  is  expected  from 
Born  approximation  results.  Jaecks,  et  al.(Ref.  12)  have  studied  the 
capture  into  the  metastable  2s  state  of  hydrogen  for  protons  incident  on 
the  rare  gases.  Considerable  work  on  capture  into  excited  projectile 
states  was  reported  at  the  1967  Leningrad  conference  (Ref.  13),  however, 
the  choice  of  projectile  still  appears  to  be  protons. 

The  two-state  ns  model  discussed  in  reference  1  and  in  Section  II-l  is 
employed  to  estimate  the  cross  section  for  capture  into  excited  states 
for  several  systems.  The  final  projectile  states  considered  are  the 
lowest  few  s  states  in  lithium.  Is -model  wavefunctlons  are  used  for  the 
target  systems  of  atomic  nitrogen  and  neon.  The  ns  states  of  the  projec¬ 
tile  are  represented  by  the  corresponding  ns  hydrogenlc  wavefunction  with 
the  effective  charge  adjusted  to  reproduce  the  experimental  ionization 
energy  of  the  given  projectile  state. 

The  results  of  the  explicit  calculations  for  lithium  ions  in  atomic 
oxygen  and  neon  target  gases  are  shown  in  Figures  1  and  2.  Attempts  were 
made  to  extend  the  calculations  to  the  capture  into  s  states  with  higher 
values  of  the  quantum  number  n,  and  for  capture  into  n  -  4  states  of 
sodium  (results  for  the  capture  into  the  n  -  3  ground  state  of  sodium  are 
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Figure 


VELOCITY  (atomic  units) 

Charge  Exchange  Croee  Section!  for  Capture  into  the 
and  n-3  State*  of  Lithium  from  Neon. 
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n-2 
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presented  in  Figures  3  and  4).  The  attempts  have  been  unsuccessful  since 
the  number  of  nodes  In  an  s  wave function  Increases  with  n,  which  Introduces 
enough  additional  oscillations  into  the  required  matrix  elements  so  that 
the  present  numerical  procedures  are  not  adequate.  It  may  be  possible  to 
calculate  capture  cross  sections  Into  excited  states  with  larger  n  than  In 
the  present  case  by  using  the  existing  programs  and  carfully  adjusting 
those  Input  parameters  that  control  the  Integration  step  length  and  the 
quadrature  formulae  for  the  matrix  elements.  An  alternate  approach  la  to 
represent  the  higher  excited  states  by  ns  model  wavefunctlons  with  lower 
values  of  n  but  with  the  effective  charge  selected  to  reproduce  the  exper¬ 
imental  value  of  the  excited  state  energy.  Such  a  procedure  can  be  ex¬ 
pected  to  yield  reliable  results  only  if  the  dominant  part  of  the  capture 
cross  section  comes  from  those  values  of  the  impact  parameter  for  which 
the  Incorrect  representation  of  the  inner  inodes  of  the  model  wavefunctlon 
are  not  important.  While  such  objections  also  can  be  made  about  the  use 
of  the  Is -model  wavefunctlons  for  the  target  and  projectile  systems  in 
the  calculation  of  the  cross  section  for  capture  into  the  ground  state,  it 
must  be  remembered  that  the  mean  radius  of  an  excited  state  orbital  la  much 
larger  than  the  mean  radius  of  the  corresponding  ground-s.tate  orbital.  From 
the  above  discussion,  it  is  clear  that  the  accurate  calculation  of  cross 
sections  is  highly  excited  states  requires  more  careful  study. 


According  to  the  Oppenhelmer  rule,  the  ratio  of  the  cross  section  for 
capture  at  high  energies  into  the  3s  state  of  lithium  to  the  cross  section 
for  capture  into  the  2s  state  should  be  (2/3)3  _  0.2963.  At  the  Impact 
velocity  of  3.0  atomic  units,  the  calculated  ratio  is  0.27  for  collisions 
with  atomic  oxygen  and  0.26  for  collisions  with  neon.  The  model  calcula¬ 
tions  show  that  the  Impact  velocity  at  which  the  cross  section  obtains  its 
maximum  value  increases  somewhat  for  the  excited  Btates  and  that  the  maxi¬ 
mum  cross  sections  depart  significantly  from  the  ratio  of  the  cubes  of  the 
quantum  numbers  n.  For  the  lithium-oxygen  cross  sections,  the  ratio  of 
the  maximum  cross  section  for  capture  into  the  ground  2s  state  to  the 
maximum  cross  section  for  capture  into  the  3s  state  is  about  4.7,  and 
the  corresponding  ratio  for  the  neon  target  gas  is  6.9. 
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Figure  31  charge  Exchange  Cross  Section  for  Capture  into  the  n-3 
State  of  Sodiua  froa  Ataalc  Oxygen. 
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SECTION  (10 


SECTION  III 


QUANTUM  THEORY  OF  ELECTRON  LOSS  IN  HEAVY -BODY  COLLISIONS 


1.  Introduction 


The  collision  processes  of  Interest  sre  of  the  type 

A  +  B-A+  +  B+  e  (28) 


where  A  and  B  are  neutral  or  charged  atoms  and  molecules.  System  A,  which 
looses  an  electron  during  the  collision, is  termed  the  projectile,  and 
system  B  Is  called  the  target  system.  It  Is  assumed  that  either  the  tar¬ 
get  or  projectile  system  Is  electrically  neutral  in  the  Initial  channel 
so  that  the  problem  la  not  complicated  by  a  long-range  Coulomb  Interaction 
Asong  the  reactants  In  the  initial  channel. 

Ionisation  processes  during  very  low-velocity  heavy-body  collisions 
have  received  some  attention  In  the  past  few  years.  The  low-energy  theo¬ 
retical  treatments  have  Involved  a  model  based  on  a  statistical  treatment 
of  the  "evaporation"  of  electrons  from  an  excited  system  created  by  the 
annual  interaction  of  the  target  and  projectile  charge  distributions,  or 
are  based  on  the  consideration  of  auto-lonlsation  or  Auger  processes. 

Same  relevant  papers  on  low-energy  Ionisation  processes  Include  Russek  . 

(Ref.  14),  Fano  and  Llchten  (Ref.  IS),  Kessel  and  Everhart  (Ref.  16),  and 
Demkov  and  Komarov  (Ref.  17),  while  Smirnov  and  Flrsov  (Ref.  18)  have  con¬ 
sidered  negative  lon-atom  systems  at  low  velocities.  The  details  of  the 
low-velocity  theories  will  not  be  reviewed  because  it  Is  unlikely  that  they 
can  be  meaningfully  extended  to  treat  moderate  and  high-energy  collisions. 

For  high-energy  heavy-body  collisions  Involving  Ionisation,  the  Belfast 
group  (Ref.  19  through  21)  have  performed  extensive  calculations,  using  the 
Born  approximation  for  light  ions  and  atoms.  Although  the  Born  approximation 
results  show  good  agreement  with  the  experimental  results  at  high  energies, 
complete  Born  approximation  calculations  for  many-e lectron  projectile  and 
target  systems  become  Involved  and  have  not  been  attempted.  Dmitriev  and 
Nikolaev  (Ref.  2)  have  proposed  a  model  to  describe  stripping  in  high-energy 
collisions  based  on  a  free-colllslon  approximation  in  which  the  Born  approx¬ 
imation  was  used  for  the  electron-target  system  scattering.  They  have  per¬ 
formed  explicit  calculations  for  light  target  and  projectile  systems  and 
find  that  the  model  yields  results  Identical  to  the  Born  approximation  at 
high  energies.  Calculations  based  on  the  free-colllslon  approximation  can 
be  extended  readily  to  many-electron  target  and  projectile  systems,  with 
little  computational  difficulties. 
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2 .  Theory 


Lee  k_  be  Che  magnitude  of  Che  Initial  momentum  rector  kQ  of  the  elec¬ 
tron  Incident  on  a  apherlcally  ayrmaetrlc  atom  with  Z  electrons.  Let  Ji,  be 
the  final  ecattered  electron  momentum,  and  K  be  the  magnitude  of  the  mo¬ 
mentum  tranafer  vector  K  ■  kj  -  .  For  elastic  acatterlng, 


e  +  B  -  e  +  B  (29) 

the  croaa  aectlon  In  momentum  apace  from  the  Born  approximation  (Ref.  22)  la 

|Z  "  F(K)|2dK  (30) 

ko  K 

where  F(K)  la  the  elaatlc  form  factor  given  by  the  ground  atate  expectation 
value 

F(K) 

>1 

the  a  urn  being  over  the  target  electrona.  For  aufflclently  large  K,  uae 
can  be  made  of  the  closure  properties  of  the  target  eigenfunctions  (Ref. 22) 
and  the  sum  of  the  cross  sections  of  all  contributing  Inelastic  processes 
of  the  type 

e  +  B  -  e  +  B(n)  (32) 

where  n  Indicates  excited  target  states.  Including  the  continuum,  can  be 
written  In  the  Born  approximation  as 


£ 


<exp  (1  K  •  £j)> 


(31) 


^In  -  7^3  2  '‘in^ 


(33) 


where  Sln(K)  la  the  Incoherent  scattering  factor  given  by 


sm<K> 


i\t< 


<exp[l  K  (rj  -  rk)]>  -  |F(K)| 


(34) 
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The  evaluation  of  Sln(K)  la  complicated  by  the  occurrence  of  the  expecta¬ 
tion  value  of  a  two-electron  operator.  In  order  to  use  the  closure  re¬ 
lationship  to  obtain  Equation  (33),  the  momentum  transfer  must  be  much 
larger  than  the  momentum  required  to  supply  the  excitation  energy  for 
all  significantly  contributing  inelastic  processes.  If  the  important 
momentum  transfer  values  are  not  large  enough  to  justify  using  the 
closure  properties,  it  would  be  necessary  to  calculate  the  differential 
scattering  cross  section  for  all  of  the  many  significant  inelastic  ex¬ 
citation  and  ionization  processes.  Such  calculations,  on  a  large  scale 
for  many-electron  targets,  are  impractical  because  of  the  need  for  accu¬ 
rate  excited  and  continuum  wavefunctions  for  the  many-electron  target 
system. 

In  the  free-colllslon  approximation,  the  projectile  system  is  assumed 
to  be  composed  of  n  electrons,  each  with  a  binding  energy  1^  required  to 
remove  the  i-th  electron,  and  each  moving  with  an  incident  momentum  kQ, 
corresponding  to  the  relative  velocity  between  the  two  nuclei.  These 
incident  electrons  are  assumed  to  suffer  elastic  and  inelastic  collisions 
as  free  electrons  with  the  target  system.  For  the  elastic  scattering, 
the  electron  la  assumed  to  be  removed  from  the  projectile  system  if  the 
momentum  transfer  it  acquires  in  the  collision  is  sufficient  to  increase 
its  energy  by  an  amount  greater  than  its  binding  energy  Ij.  For  the 
inelastic  scattering,  the  electron  is  assumed  to  be  removed  if  it  acquires 
enough  energy  through  the  momentum  transfer  to  exceed  both  the  binding 
energy  of  the  electron  to  the  projectile  and  the  excitation  energy  re¬ 
quired  to  excite  the  target.  The  total  stripping  cross  section  la  given 
by 


where  the  sum  is  taken  over  the  projectile  electrons  and  the  limits  of 
integration  are  determined  by  the  kinematics.  The  lower  limit  for  the 
elastic  contribution,  k^  is  given  by 

ku  -  -J2  It  (36) 

and  the  upper  limit  is  given  by 

k21  -  2  ko  (37) 


The  lower  limit  for  the  inelastic  contribution  is 
quantities 


and 


the  greater  of  the  two 
-  2AE/k02)*] 


16 


where  the  first  quantity  represents  the  minimum  momentum  transfer  required 
to  Increase  the  energy  of  the  1-th  target  electron  by  Iit  and  the  second 
quantity  Is  the  minimum  momentum  trasfer  needed  to  trasfer  an  excitation 
energy  AE  to  Che  target  system.  Dmitriev  and  Nikolaev  (Ref.  2)  suggest 
using  I0>  the  first  ionization  potential  of  the  target  system  as  AE.  At 
high  values  of  the  Incident  momentum  k0,  l.e.,  kQ2  »  AE,  the  choice  of 
AE  is  unimportant  because  the  first  quantity  Is  Independent  of  k„  and  the 
second  decreases  as _AE/kQ  so  that  for  large  k9>  the  lower  limit  Is  deter¬ 
mined  entirely  by  V2  I i *  except  possibly  for  Inner  shell  electrons, which 
make  little  contribution  to  the  sum  in  Equation  (35).  For  processes  In 
which  AE  >  l£,  for  example  the  stripping  of  the  loosely  bound  electron 
from  a  negative  ion,  the  value  of  AE  is  important  because  the  lower  limit 
is  determined  by  the  expression  containing  AE  for  a  considerable  range  of 
initial  momentum  kg.  For  this  energy  range, the  use  of  Equation  (33)  Is 
not  valid  because  the  conditions  necessary  for  the  use  of  the  closure 
properties  Is  not  justified.  It  has  been  assumed  that  AE  Is  given  by  the 
logarithmic  mean  energy , which  enters  In  the  Bethe  (Ref.  23)  theory  for 
the  stopping  of  fast  heavy  particles  In  matter,  so  that 


in  AE 


S  foj  *"(EoJ> 


(38) 


where  fQ.  Is  the  oscillator  strength  for  excitation  of  the  target  system  from 
the  ground  state  to  the  level  J,  E0j  is  the  corresponding  excitation  energy, 
and  the  summations  in  Equation  (38)  include  integrations  over  the  continuum. 

The  use  of  Equation  (38)  to  determine  Che  lower  limit  of  the  inelastic  scat¬ 
tering  integration  can  be  considered  as  using  an  effective  upper  limit  for 
AE.  The  upper  limit  for  the  Inelastic  scattering  contribution  must  also 
be  selected  with  care.  From  purely  klnematlcal  arguments,  k^  for  the 
inelastic  electron-atom  scattering  would  be  about  2  kg.  However,  Equation  (33) 
fails  at  very  high  energies  and  the  maximum  momentum  transfer  Involved  Is  the 
maximum  momentum  transfer  In  electron-electron  scattering  processes  which  Is 
kQ.  According  to  Dmitriev  and  Nikolaev  (Ref.  2),  the  following  is  used 
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■  k 


(39) 


Both  the  form  factor  F(K)  and  the  Incoherent  scattering  factor  Sin(K) 
are  Important  in  the  field  of  X-ray  physics,  the  former  in  elastic  X-ray 
scattering  and  the  latter  in  Compton  scattering  of  X-rays,  so  that  con¬ 
siderable  work  on  the  calculation  of  the  scattering  factors  Is  available 
from  the  literature.  Kim  and  Inokutl  fRef.  24)  have  calucated  highly  ac¬ 
curate  values  of  F(K)  and  S^n(K)  for  helium  by  using  the  accurate  correlated 
wavefunction  of  Hart  and  Herzberg  (Ref.  25).  For  systems  with  more  elec¬ 
trons,  Hartree-Fock  calculations  have  been  performed  for  the  scattering 
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factors,  •  recent  tabulation  being  reference  26.  Many  of  the  values  of  the 
scattering  factors  given  In  reference  26  are  taken  from  the  work  of  Freeman 
(Ref.  27),  who  considers  systems  with  other  than  S  ground-state  syimnetry. 
Cromer  and  Mann  (Ref.  28)  have  recently  published  Incoherent  scattering 
factor  calculations  for  most  closed  shell  systems,  using  Hartree-Fock 
wavefunctlons .  For  the  target  systems  of  Interest  In  this  study.  It  Is 
not  necessary  to  calculate  scattering  factors  because  the  values  have  been 
tabulated  In  the  literature.  For  nonspherlcal  target  systems,  the  spheri¬ 
cal  target  systems,  the  spherelcal  average  scattering  factors  have  been 
employed  to  be  consistent  with  the  neglect  of  the  symmetry  properties  of  the 
projectile  system.  For  all  atomic  systems,  the  elastic  scattering  factor 
Is  a  monotonlca 1 ly  decreasing  function  of  K,  having  the  value  Z  for  K  »  0 
and  decreasing  to  zero  for  large  K.  The  Incoherent  scattering  factor  is 
a  monotonies lly  Increasing  function  of  K,  assuming  the  value  0  at  K  «  0, 
and  Increasing  to  Z  for  large  values  of  K.  The  Integrands  for  both  the 
elastic  and  Inelastic  contributions  to  Equation  (35)  are  smooth  functions 
of  K  so  that  the  Integrals  are  performed  simply  by  standard  numerical 
methods.  The  Integrands  for  helium, using  the  calculations  of  Kim  and 
Ionkutl  (Ref.  24),  are  shown  in  Figure  5.  From  the  behavior  of  the  In¬ 
elastic  lnteg^«nd  for  low  values  of  momentum  transfer,  it  can  be  seen 
that  the  Inelastic  contribution  Is  more  sensitive  to  the  value  of  the 
lower  limit  of  Integration  that:  the  elastic  contribution. 

Both  Integrands  decrease  rapidly  for  large  K,  so  that  the  large  kQ, 
the  upper  limits  may  be  replaced  by  »  and  the  high  -energy  behavior  of 
the  stripping  cross  section  becomes 


a/k02  (40) 


|Z  -  F(K)|2  +  Z  Sln(K)  (41) 


The  E_1  high-energy  dependence  of  the  stripping  cross  section  is  consistent 
with  experimental  data  (see  Section  1V-2). 

The  information  needed  for  the  projectile  system  Is  the  number  of  elec¬ 
trons  In  each  shell  and  the  one-electron  energies  required  to  remove  the 
electron  from  the  individual  shells.  For  neutral  and  positive  Ion  projec¬ 
tiles,  the  energy  needed  to  remove  the  outermost  electron  Is  the  Ionization 
potentials , which  are  taken  from  the  experimental  tabulations  of  Moore  (Ref. 
29).  For  the  inner  Bhells,  the  estimates  of  Sister  (Ref.  30)  are  used. 
Results  of  detailed  calculations  and  comparison  with  experimental  results 
are  presented  In  Section  IV-2 . 
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SECTION  IV 
RESULTS 


1.  Choree  Exchange 

The  results  of  using  Che  baalc  la  model  are  presented  In  Tables  II 
through  VIII  and  graphed  In  Figurea  6  through  46  where  the  experimental 
reaulta  of  Layton  and  Fite  (Ref.  31)  are  alao  shown.  The  calculated 
croaa  sections,  as  might  be  expected  from  such  a  simple  model,  share  the 
saaM  general  shape  although  the  predicted  maximum  cross  sections  vary 
through  several  orders  of  magnitude  for  different  reactions.  In  the 
case  of  oxilecular  O2  and  N2  targets,  these  have  been  treated  as  single 
systeau  with  approximately  the  same  ionisation  potentials  as  the  atomic 
aysteaw,  thus  a  factor  of  two  has  not  been  applied  to  the  calculated 
cross  sections  In  these  cases. 

The  la  model  cross  sections  are  usually  amaller  than  the  experlaientally 
measured  values.  Nhere  a  comparison  with  the  experimental  results  la  pos¬ 
sible,  the  la  modal  yields  a  sharper  amxlmum  than  the  former,  from  there 
falling  off  much  too  rapidly  as  the  energy  diminishes  and  soaiewhat  more 
rapidly  In  aoam  caaaa  as  the  energy  increases.  Use  of  the  n»  model  may 
give  a  leas  rapid  fall-off  at  the  hlgh-veloclty  and, as  may  be  seen  from 
Figure  47.  At  a  sufficiently  high  energy, the  exchange  correction  will 
certainly  operate  In  the  favorable  sense. 
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CROSS  SECTIONS  CT,_  (UNITS  OF  10'17  cm2)  CALCULATED  FOR  K+  PROJECTILES  (is  MODEL) 
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The  number  In  parentheses  is  the  power  of  10  by  which  the  entries  are  to  be  multiplied. 


CROSS  SECTIONS  a,,  (OMITS  OF  10  cm)  CALCULATED  FOR  I  PROJECTILES  (Is  MODEL) 
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The  number  In  parentheses  is  the  power  of  10  by  which  the  entries  are  to  be  multiplied. 


CROSS  SECT  IOCS  o’  (UNITS  OF  10"  cm  )  CALCULATED  FOR  Cr  PROJECTILE  (is  MODEL) 
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The  nuaber  in  parentheses  is  the  power  of  10  by  which  the  entries  are  to  be  aultiplied. 
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CROSS  SECTIONS  a,. (UNITS  OF  10  cm  )  CALCULATED  FOR  N*  PROJECTILE  (Is  MODEL) 
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CROSS  SECTIONS  7  (UNITS  OF  10  ca  )  CALCULATED  FOR  Rb  PROJECTILE  (Is  MODEL) 
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CROSS  SECTIONS  a..  (UNITS  OF  10  cm)  CALCULATED  FOR  Cs  PROJECTILE  (Is  MODEL) 
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WOSS  SECTION 


Figure  8.  Charge  Exchange  Cron  Section  foi  Capture  Into  the  Ground  State 
for  Singly  Charged  Fotaeslum  In  Neon. 
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Figure  9.  Charge  Exchange  Croce  Section  for  Capture  into  the  Ground  State 
for  81nglp  Charged  rotaeelun  in  Argon.  The  Experimental  Raeults 
ere  from  Inference  31. 
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Figwra  10.  Chart*  tech— t*  Cross  Section  for  Cap  Curs  into  eh*  Ground  Scat* 

for  Stilly  Charged  Potaaslua  In  Atosie  Mltrogsn.  Ths  tepsrl— ntal 
tesults  srs  froa  Ksfsrsnea  31. 
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F lgure  12.  Charg*  Exchange  Croaa  Sac t Ion  for  Capture  InCo  tha  Ground  State 
for  Singly  Charged  Iodine  In  Atonic  Hydrogen. 


Figure  15.  Charge  Exchange  Cross  Section  for  Capture  into  the  Ground,  State 
for  Singly  Charged  Iodine  in  Argon.  The  Experimental  Results 
are  from  Reference  31. 
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Figure  16.  Charge  Exchange  Cross  Section  for  Capture  Into  the  Ground  State 
for  Singly  Charged  Iodine  In  Atonic  Nitrogen.  The  Experimental 
Results  are  from  Reference  31. 
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Figure  17.  Charge  Exchange  Croas  Section  for  Capture  Into  the  Ground  State 
for  Singly  Charged  Iodine  In  Atomic  Oxygen.  The  Experimental 
Reaulte  are  from  Reference  31. 
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Figure  22.  Charge  Exchange  Croas  Section  for  Capture  Into  the  Ground  State 
for  Singly  Charged  Chromium  in  Atomic  Nitrogen. 
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Figure  23.  Charge  Exchange  Croat  Section  for  Capture  into  the  Ground  State 
for  Singly  Charged  Chroalum  in  Itoelc  Oxygen. 
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Figure  26.  Charge  Exchange  Cross  Section  for  Capture  Into  the  Ground  State 
for  Singly  Charged  Lithium  In  Neon. 
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Figure  27.  Charge  Exchange  Cross  Section  for  Capture  into  the  Ground  State 
for  Singly  Charged  Lithium  in  Argon. 
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Figure  28.  Charge  Exchange  Cross  Section  for  Capture  into  the  Ground  State 
for  Singly  Charged  Lithium  in  Atomic  Nitrogen. 
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Figure  30.  Charge  Exchange  Croas  Section  for  Capture  Into  the  Ground  State 
for  Singly  Charged  Sodium  In  Atomic  Hydrogen. 
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Figure  37.  Charge  Exchange  Cross  Section  for  Capture  into  the  Ground  State 
for  Singly  Charged  Rubidium  in  Helium. 
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Figure  40.  Charge  Exchange  Cross  Section  for  Capture  into  the  Ground  State 
for  Singly  Charged  Rubidium  in  Atomic  Oxygen. 
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Charge  Exchange  Cross  Section  for  Capture  into  the  Ground  S 
for  Singly  Charged  Cesium  in  Argon. 
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SECTION 


2 .  Electron  Stripping  Cross  Sections 


The  theory  developed  In  Section  III  has  been  employed  to  calculate 
the  electron  stripping  cross  sections  for  a  variety  of  ions  in  several 
atomic  gases.  The  projectile  systems  considered  are  H,  H" ,  He  He+,  Li, 

Li",  N,  0,  Ne ,  F,  Na ,  Mg,  A£,  Fe,  Mn,  and  Cr .  The  target  gases  con¬ 
sidered  are  He,  N,  0,  Ne ,  and  Ar.  A  wide  range  of  impact  energies  are 
considered,  however,  as  discussed  in  Section  III,  the  theory  is  limited 
to  high  energy  impacts.  In  this  section,  the  results  of  these  calcula¬ 
tions  are  presented  and  the  theoretical  predictions  are  compared  with 
experimental  data. 

In  Figure  48, the  stripping  cross  section  is  presented  for  the  few- 
electron  system  process 

H  +  He  -  H+  +  e  +  ^  He  (42) 

where  by  ^  He  the  elastic  and  inelastic  processes  is  indicated  including 

ionization  of  the  helium  atom.  The  experimental  data  of  Allison  (Ref.  32), 
Berkner ,  et  al. (Ref.  33),  and  Smythe  and  Toevs  (Ref.  34)  are  shown  for 
comparison.  At  high  energies,  the  cross  section  is  given  by 


(43) 


-11  2 

with  a  -  1.595  x  10  cm  eV  and  E  is  the  projectile  energy  in  the  labora¬ 
tory  frame  in  electron  volts.  The  agreement  between  theory  and  experiment 
is  good  over  the  energy  range  of  100  keV  to  100  MeV,  and  is  remarkable  at 
the  10  MeV  and  14.6  MeV  data  points.  Figures  49  and  50  give  the  results 
for  hydrogen  atoms  incident  on  the  heavier  system  of  atomic  nitrogen  and 
argon.  For  these  target  systems,  the  cross  section  at  high  energies  can 
again  be  represented  by  Equation  (43)  with  a  assuming  the  values  4.237x10” 
cm2  ey  and  4.706  x  10"^®  cm^  eV  for  nitrogen  and  argon,  respectively.  The 
agreement  between  theory  and  experiment  is  not  as  good  as  for  the  H-He  case. 
At  100  keV  the  theory  overestimates  the  cross  section  by  a  factor  of  about 
3  for  the  nitrogen  target,  where  the  per  atom  result  for  molecular  nitrogen 
is  used  for  the  experimental  data,  and  a  factor  of  about  7  for  the  argon 
target.  At  higher  energies,  where  the  theory  is  valid,  the  agreement  is 
better;  the  theory  giving  results  about  10  and  20  percent  higher  than  the 
experimental  results  for  the  nitrogen  and  argon  targets,  respectively. 

Kim  and  Inokuti  (Ref.  24)  have  used  Hartree-Fock  wavefunction  as  well 
as  accurate  correlated  wavefunctions  in  their  calculations  of  the  atomic 
form  factor  and  the  incoherent  scattering  function.  These  results  give 
an  atomic  form  factor  using  the  correlated  wavefunction  that  is  larger 
for  large  K  values  and  smaller  at  small  K  values  than  the  corresponding 
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Figure  48.  Electron  Stripping  Cross  Section  for  Hydrogen  Atoms  Incident  on 
Atomic  Hydrogen.  The  Experimental  Results  are:  O  .Ref.  32; 

Ref.  33  ;A,  Ref.  34. 
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CROSS  SECTION 


Figure  49-  Electron  Stripping  Cross  Section  for  Hydrogen  Atoms  Incident  on 

Atomic  Nitrogen.  The  Experimental  Results  are:  Q  ,  Ref.  32 ;  (  |  , 
Ref.  33;  A  ,  Ref.  34. 


72 


Figure  50.  Electron  Stripping  Cross  Section  for  Hydrogen  Atoms  Incident  on 
Argon.  The  Experimental  Results  are:  O.  Ref.  32;  Ref.  33; 
A  ■  Ref.  34. 
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Hartree-Fock  results.  Hie  more  accurate  calculations  of  the  incoherent 
scattering  factors  using  the  correlated  wavefunctions  is  smaller  than  the 
corresponding  Hartree-Fock  results  for  all  values  of  the  momentum  trans¬ 
fer.  It  is  expected  that  this  behavior  of  the  scattering  factors  is  not 
limited  to  the  helium  system,  but  is  true  for  larger  systems.  As  the 
number  of  electrons  in  the  system  increases,  the  accuracy  of  the  Hartree- 
Fock  results  decreases,  especially  for  the  Incoherent  scattering  factor, 
which  is  a  two-electron  operator  expectation  value.  Because  the  Hartree- 
Fock  Incoherent  scattering  factor  is  larger  at  all  K  than  the  accurate 
Incoherent  scattering  factor,  the  integrand  for  the  inelastic  scattering 
contribution  is  overestimated  by  the  use  of  Hartree-Fock  incoherent 
scattering  factors  at  all  K,  so  that  the  contribution  of  inelastic  pro¬ 
cesses  to  the  high-energy  stripping  cross  section  is  overestimated  by 
the  Hartree-Fock  theory.  For  the  elastic  contribution,  similar  arguments 
show  that  for  large  K  the  use  of  Hartree-Fock  form  factors  overestimates 
the  integrand,  but  at  small  K  it  underestimates  the  integrand.  The  net 
effect  on  the  elastic  contribution  can  vary  from  system  to  system,  but  at 
high  energies  for  projectile  systems  with  large  ionization  potentials, 
where  the  lower  limit  of  the  elastic  integration  is  large,  the  use  of 
Hartree-Fock  form  factors  will  also  lead  to  an  overestimation  of  the 
stripping  cross  section.  The  extent  to  which  the  use  of  Hartree-Fock 
scattering  factors  contributes  to  t he  disagreement  between  theory  and 
experiment  for  the  high-energy  H-N  and  H-A  r  results  cannot  be  resolved 
because  it  is  dependent  on  the  availability  of  more  accurate  nitrogen  and 
argon  wavefunctions ,  which  must  await  considerable  theoretical  and  com¬ 
putational  advances.  However,  some  of  the  discrepancy  is  caused  by  the  use 
of  Hartree-Fock  wavefunctions  in  the  scattering  factor  calculations. 

Figure  51  compares  the  theoretical  cross  section  for  collisional  de¬ 
tachment  of  the  electron  from  the  negative  hydrogen  ion  upon  impact  with 
helium  with  the  experimental  data  from  Allison  (Ref.  32).  For  these 
impact  energies,  the  theory  is  not  expected  to  apply  so  that  the  agree¬ 
ment  to  within  a  factor  of  about  two  may  be  accidental.  The  process  is 
so  efficient  that  the  cross  section  has  not  been  measured  accurately  at 
very  high  energies  because  the  equilibrium  fraction  of  negative  ions  in 
the  beam  becomes  very  small.  It  is  interesting  to  note  that  here  the 
theoretical  results  lie  below  the  experimental  results. 

Figure  52  shows  the  electron  stripping  cross  section  for  the  many- 
electron  projectile  system  of  sodium  incident  on  helium  gas.  The  con¬ 
tributions  of  the  various  shells  to  the  total  cross  section  is  shown. 

At  high  energies,  it  can  be  seen  that  the  2p  shell  contribution  is 
significantly  larger  than  the  3s  shell  contribution  although  the  energy 
required  to  reirove  the  2p  electrons  is  larger  than  the  energy  required 
to  remove  the  3s  electrons.  At  high  energies,  the  ratio  of  the  2p  con¬ 
tribution  to  the  3s  contribution  is  roughly  3:1,  while  the  ratio  of  the 
number  of  2p  electrons  to  3s  electrons  is  6,  so  that  per  electron,  it 
is  more  efficient  to  remove  the  outer  electron.  The  contribution  of  the 
2s  shell  is  negligible,  as  is  the  Is  contribution.  The  abrupt  changes 
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Figure  51.  Collisional  Detachment  Cross  Section  for  Negative  Hydrogen  Ions 
Incident  on  Helium.  The  Experimental  Results  are  from  Ref.  32. 
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Figure  52.  Electron  Stripping  Cross  Section  for  Sodium  Ions  on  Helium  Showing 
the  Total  Cross  Section  and  the  Contributions  of  the  Different 
Shells . 
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of  the  cross  sections  shown  are  not  real  and  are  only  an  artifact  of  the 
theory.  At  \e  cusps,  the  lower  limit  of  the  integration  for  the  inelastic 
contribution  changes  from  being  dependent  on  kD  and  AE,  to  being  a  constant 
dependent  only  on  the  one  electron  energies  of  the  electrons.  However, 
the  magnitude  of  the  variation  in  the  vicinity  of  the  cusps  show  that  the 
inelastic  scattering  contribution  is  large  and  is  a  sensitive  function  of 
the  lower  limit  of  the  integration  as  discussed  previously  in  Section  III. 

Some  calculations  of  the  electron  stripping  cross  section  for  a  more 
highly  stripped  projectile  system  have  been  performed  involving  an  argon 
target  and  Cr ,  Cr+,  Cr"1"*,  Cr+3,  Cr-^,  and  Cr^projectiles .  The  results 
are  shown  in  Figure  53.  At  low  energies,  the  cross  section  decreases 
rapidly  with  the  stage  of  ionization  of  the  projectile,  but  for  higher 
energies  the  decrease  in  cross  section  is  not  nearly  so  rapid,  and  the 
highly  stripped  ions  still  have  a  large  cross  section  for  losing  at. 
additional  electron.  At  these  high  energies, the  equilibrium  fraction  of 
highly  stripped  ions  in  a  beam  is  large.  For  the  more  highly  stripped 
ions  included  in  Figure  53,  accurate  values  for  the  energy  to  remove  one 
electron  for  the  shells  just  below  the  valance  shell  are  not  available, 
so  that  estimates  were  employed.  Accurate  values  of  the  valence  electron 
ionization  potentials  (Ref.  29)  are  available,  however,  and  in  Figure  54 
only  the  3d  valence  shell  contribution  to  the  stripping  cross  section  is 
shown.  These  results  can  be  considered  as  lower  limits  to  the  stripping 
cross  section  in  the  free-electron  scattering  approximation.  These 
values  show  a  more  rapid  decrease  of  the  stripping  cross  section  as  the 
stage  of  ionization  of  the  projectile  increases  than  the  corresponding 
results  when  inner  shell  contributions  are  included. 

The  results  of  the  calculations  for  che  other  systems .which  are  pre¬ 
sented  in  Figures  55  through  64, will  not  be  discussed  in  detail.  Because 
of  the  limitations  of  the  theory,  only  the  results  at  very  high  energy 
are  expected  to  be  qualitatively  reliable. 
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Figure  55.  Electron  Stripping  Cross  Sections  for  Various  Projectile  Systems  in  Helium  Gas. 
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Figure  56.  Electron  Stripping  Cross  Sections  for  Various  Projectile  Systems  in  Helium  Gas. 
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Figure  57.  Electron  Stripping  Cross  Sections  for  Various  Projectile  Systems  in  Atomic 
Nitrogen  Gas. 
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Figure  59.  Electron  Stripping  Cross  Sections  for  Various  Projectile  Systems  in  Atomic 
Oxygen  Gas. 
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Figure  61.  Electron  Stripping  Cross  Sections  for  Various  Projectile  Systems  in  Neon 
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for  Various  Projectile  Systems  in  Neon. 
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Figure  63.  Electron  Stripping  Cross  Sections  for  Various  Projectile  Systems  In  Argon 
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SECTION  V 


ENERGY  LOSS  OF  FAST  HEAVY  IONS  IN  GASES 


The  Bethe  (Ref.  23)  theory  of  the  stopping  of  fast  charged  heavy 
particles  In  neutral  gases  has  been  successful  in  the  interpretation  of 
the  high  Impact  energy  data.  The  Bethe  theory  is  based  on  the  Born  approx¬ 
imation  with  certain  additional  assumptions  and  ignores  the  possibility 
of  excitation,  electron  capture,  and  electron  loss  by  the  projectile  system. 
The  neglect  of  inelastic  processes  for  the  projectile  and  the  use  of  the 
Born  approximation  cause  the  Bethe  theory  to  fall  at  low  energies. 

Dalgarno  and  Griffing  (Ref.  35)  have  performed  a  detailed  study  of  the 
rate  of  energy  loss  of  protons  in  atomic  hydrogen.  They  list  the  possible 
inelastic  processes  and  obtain  expressions  for  the  contribution  to  the 
total  stopping  power  for  each  process.  The  total  stopping  power  is  ex¬ 
pressed  in  terms  of  the  equilibrium  fractional  charge  content  of  the  beam 
and  the  energy  loss  cross  sections  for  the  various  processes.  Explicit 
results  using  the  Born  approximation  cross  sections  for  the  relevant  pro¬ 
cesses  are  given. 

The  formulation  of  the  energy  loss  theory  of  Dalgarno  and  Griffing 
can  be  extended  straightforwardly  to  beam  and  target  systems  with  more 
electrons.  If  it  is  assumed  that  the  stopping  gas  consists  entirely  of 
an  atomic  gas  Y  in  its  ground  state,  the  inelastic  processes  that  must  be 
considered  in  detailed  energy  loss  calculations  are 

*W+  nf+  mt+ 

X  (j)  +  Y  -  X  (k)  +  Y  (£)  +  (mf  +  nf  -  n^e  (44) 

where  n^,  nf,  and  mf  are  positive  or  negative  integers  (if  the  negative 
ions  exist)  or  zero,  and  j,  k,  and  l  represent  the  quantum  numbers  neces¬ 
sary  to  specify  the  possible  states  of  the  systems.  If  a  proton  beam 
Incident  on  atomic  hydrogen  is  considered,  the  range  of  values  that  nf, 
n j ,  and  mf  can  assume  is  small  because 

-1  <  nf  <  +1 

-i  <  nf  <  +i  (45) 

-1  <  mf  <  +1 
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and  only  one  state  exists  for  the  positive  and  negative  ion.  For  many- 
electron  target  and  projectile  systems  at  high  energies,  more  highly 
stripped  projectile  ions  can  become  an  important  fraction  of  the  equili¬ 
brium  beam,  and  the  number  of  reactions  that  must  be  considered  becomes 
large.  At  sufficiently  high  energies,  multiple  ionization  processes  may 
become  significant. 

The  electron  stripping  cross  sections  decrease  much  less  rapidly  with 
increasing  projectile  energy  than  the  charge  transfer  cross  sections.  Be¬ 
cause  of  the  rapid  decrease  of  the  charge  exchange  cross  sections  with 
increasing  energy,  at  high  energies  where  the  stripping  cross  sections  cal¬ 
culated  by  using  the  free-scattering  approximation  are  valid,  the  electron 
loss  cross  section  for  a  typical  neutral  or  once  ionized  projectile  system 
is  much  larger  than  the  corresponding  electron  capture  cross  section,  as 
indicated  by  the  results  given  in  Section  IV.  To  calculate  the  equilibrium 
charge  fractions  of  the  beam,  it  would  be  necessary  to  calculate  charge 
exchange  cross  sections  for  highly  stripped  projectile  ions  incident  on 
the  neutral  target  gas.  Because  such  cross  sections  are  not  presently 
available,  mer-viingful  calculations  for  the  rate  of  energy  loss  at  very 
high  energies  cannot  be  performed.  Similarly,  meaningful  calculation  for 
the  rate  of  energy  loss  at  lower  energies  must  await  the  accurate  deter¬ 
mination  of  ionization  and  stripping  cross  sections  in  the  desired  energy 
range.  Excitation  of  the  target  system  does  not  enter  into  the  equations 
that  determine  the  equilibrium  charge  fractions  of  the  beam,  and  excita¬ 
tion  of  the  projectile  systems  Influences  the  beam  equilibrium  only  if 
the  stopping  gas  is  sufficiently  dense  or  important  metastable  species 
formed.  However,  excitation  processes  make  significant  contributions  to 
the  energy  loss  cross  section  at  most  impact  energies,  and  must  be  known 
in  order  to  perform  meaningful  energy  loss  calculations. 

The  contribution  to  the  energy  loss  cross  section  because  of  the  charge 
transfer  process 


X+  +  Y  -  X(i)  +  Y+(m) 


(46) 


(where  X  and  Y  are  not  identical)  is 


V  '  V 1 

>  )  E  Q 

/  .  /  ■  fm  xfm 

C.  m 


where  Q is  the  cross  section,  and  the  energy  E  is  given  by 


E£m  =  T  +  E[X(£) ]  +  E[YT(m)]  -  E(XT)  -  E(Y) 


(47) 
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where  E(Z)  is  the  energy  of  the  isolated  system  Z  and  T  is  the  kinetic 
energy  acquired  by  the  captured  electron.  At  moderate  and  high  energies, 
the  kinetic  energy  term  T  constitutes  a  significant  part  of  E,  . 

iJlTl 

The  contribution  to  the  energy  loss  cross  section  from  the  stripping 
reaction 


X  +  Y  -  x"J"  +  Y  +  e  (e  ) 

m 


(48) 


Where  the  ejected  electron  has  a  kinetic  energy  e  with  respect  to  the 
parent  nucleus,  and  the  final  state  of  the  target  system  is  the  ground 
state,  is  given  by 


€  (m ) 

r  max N  ' 

E 

J  Qm(6)  (Im  +  e)d£ 

(49) 

m 

J  o 

where  Qm(€)  Is  the  cross  section  for  the  ejection  of  an  electron  with 
energy  between  e  and  e+de  and  €tnax(m)  is  determined  by  the  conservation 
of  energy  and  momentum.  Here  1^  is  the  sum  of  the  ionization  potential 
of  the  projectile  and  the  excitation  potential  of  the  final  projectile 
state  m. 

In  the  free-scattering  approximation,  only  those  excited  states  of 
the  projectile  that  can  be  represented  by  the  removal  of  a  single  valence 
or  inner  shell  electron  in  a  central  field  approximation  are  considered. 
The  energy  loss  cross  section  for  reaction  (48)  is  approximated  by 


_1 

2 


dcr 


el 


(50) 


where  the  variables  are  defined  in  Section  III.  Expression  (50)  can  be  con¬ 
sidered  a  good  approximation  to  expression  (49)  only  if  the  major  contribu¬ 
tion  to  the  Integral  in  expression  (49)  comes  from  the  part  of  energy  space 
where  e  >  I,*,  so  that  reliable  results  are  expected  only  at  high  energies 
for  systems  where  the  stripped  electron  is  loosely  bound  (1^,  small). 


The  contribution  to  the  energy  loss  cross  section  from  the  reactions 


and 


X  +  Y  -  x£  + 

e(e)  -f  Yn 

(51) 

X  +  Y  -  X+  +  Y+  + 
m  n 

e(e)  +  e(e') 

(52) 
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cannot  be  obtained  with  similar  accuracy  from  the  free-scattering  approxi¬ 
mation  because  of  the  use  of  the  closure  properties.  The  contribution  to 
the  energy  loss  cross  section  from  reaction  (51)  is 

(m,n) 

Qm,n(e>  (Im  +  En  +  e>d*  <53> 

m  n 


where  Qm  n(e)  is  the  cross  section  and  En  is  the  excitation  energy  of  the 
final  target  state.  The  contribution  to  the  energy  loss  cross  section 
from  reactions  (52 )  is 


de'  V£-e')(I»+En+£+£,)  (54) 


where  Q,,,  n(e,e')  is  the  corresponding  cross  section.  By  using  the  closure 
properties,  the  sum  of  expressions  (53)and  (54)  is  approximated  by 


i 


1 

2 


d  ct 


in 


(55) 


The  use  of  expression  (55)  is  valid  only  if  the  major  part  of  the  contribu¬ 
tion  to  the  integrals  in  expression  (53)  comes  from  the  part  of  the  energy 
space  where  e  >  IjJ,  +  En,  and  similarly  for  expression  (54)  where  the  con¬ 
dition  is  e  >  I  *  +  E  +  e'- 
m  m 

The  sum  of  expressions  (50)  and  (55)  have  been  calculated  in  units  of 
eV  cm^ ,  for  various  projectile  ions  in  argon  gas.  The  results  are  shown 
in  Figures  65  and  66.  The  approximate  energy  loss  cross  section  maxima 
occur  at  higher  Impact  energies  than  the  corresponding  e’ectron  stripping 
cross  section  maxima,  and  at  large  energies  they  show  more  gradual  decrease 
with  increasing  Impact  energy. 
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Figure  65.  Cross  Section  for  Energy  Loss  in  Electron  Stripping  Processes  for  Various  Ions  in 
Argon  Gas. 
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APPENDIX  I 


PROGRAM  XSCTN2  AND  ASSOCIATED  SECTIONS 


Several  modifications  have  been  made  to  the  program  discussed  in 
reference  1.  A  description  of  the  present  version  of  the  program  is 
included  herein.  The  basic  operation  of  the  program  remains  unchanged. 

The  modifications  affect  mainly  the  Input  and  output  and  the  option  of  a 
correction  for  the  long-range  behavior  of  the  h£i  and  hjj  matrix  elements. 
Arrangements  for  timing  the  program's  execution  have  also  been  made.  An 
uncompleted  modification  to  incorporate  the  "exchange"  correction  is  in¬ 
cluded,  marked  by  some  parameters  noted  as  "unused."  These  parameters 
should  be  provided. 

XSCTN2  computes  the  charge  exchange  cross  section  by  organizing  inte¬ 
grations  (MLTI04)  of  the  two-state  equations  along  a  series  of  trajectories 
on  the  assumption  that  the  capture  probabilities  behave  as  exp(-7p).  >(GMM) 
is  found  by  specifying  two  large  Impact  parameters  pa(RHOA)  and  pj,(RHOB)  at 
which  the  numerically  derived  capture  probabilities  Cn(pa>  £A)  and  Cn(pb) 
(CB)  are  computed  so  that 


£n  [Cn(pb)/Cn(Da)] 

^  <P."-"Pb> 

or  else  if  pa  is  set  zero,  y  is  taken  as  pj,.  y  is  assumed  to  be  the  same 
for  the  capture  probabilities  obtained  numerically  (Cn) ,  by  the  Rapp  - 
Francis  expression  (Cr)  and  by  application  of  the  model  (Cm) .  A  series 
of  pivotal  impact  parameters  is  then  generated  with  which  to  compute 
the  cross  section  quadratures  of  Gauss -Laguerre  type: 

Q  -  2n  f  dp  p  exp(-7P>  [C  exp(7p)] 

•'o 

“  2nX  wi  C^pi ^  exP(7'P1) 
i 

The  first  10  terms  of  a  16-term  quadrature  are  employed  to  compute  Q  , 

Q  ,  and  Q  corresponding  to  C  ,  C  ,  and  C  . 
xr  m  r  n  r  m 

MLTI04  performs  the  integrations  along  the  trajectory,  numerically 
(calls  INTGRT)  and  analytically,  by  evaluating  the  Rapp-Francis  expression 
and  if  MODE  =  2,  the  analytic  form  of  the  solution.  MLTI04  uses  the  as¬ 
sumed  maximum  distance  (EFFRAD)  of  interaction  to  determine  trajectory 
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length  and  also  establishes  a  "logarithmically"  spaced  set  of  reporting 
"stations"  along  the  trajectory  (clustered  more  closely  about  the  point  of 
closest  approach)  df  which  significant  quantities  in  the  calculation  are 
stored  and  from  which  are  subsequently  written.  If  MODE  =  3  these  inter¬ 
mediate  quantities  are  not  printed. 

INTGRT  performs  the  numerical  integration  along  the  trajectory  by  using 
Hamming's  method.  When  entered,  INTGRT  performs  a  step  in  the  solution  of 
the  equation  by  calling  HAMMING, 


$2.  = 


dx 


-  [c]  y 


where  [C]  y  is  provided  by  calling  RHS,  which  in  turn  calls  QCM  to  generate 
the  matrix  [C]  at  the  required  point.  INTGRT  is  used  not  only  to  integrate 
the  two-state  Equations  of  Motion,  but  also  to  integrate  two  additional 
quantities  with  respect  to  distance  along  the  trajectory,  it  adjusts  the 
Interval  of  integration.  QCM  organizes  the  subroutines  W  (through  NTRNS) 
and  QG1Q  to  perform  respectively  the  nontrans lational  and  translational 
integrals  in  the  evaluation  of  the  matrix  elements.  In  the  construction 
of  the  basic  matrix  elements  each  cross  product  arising  from  the  expressions 
for  wavefunctlons  and  potential  is  treated  separately. 

XYZ  sets  up  the  coordinate  dependent  Cll  and  C22  arrays  for  the  non- 
trans lational  integral  evaluation  performed  by  W,  X,  Y,  and  Z. 

QGLG  performs  the  quadrature  required  in  the  translational  integrals 
for  Sjj,  hji,  and  h,^  in  which  the  auxiliary  function  is  computed  by  QXF  by 
using  the  array  AXFN.  This  array  is  coded  as  follows:  the  auxiliary 
function  fni,n2  identified  by: 

J  «  r*2  +  8n^  +  1 

The  N-th  term  of  the  integrand  in  the  quadrature  with  respect  to  x  is 


and  is  give'*  by : 


C 

n 

i 


i*C 

xn  VO 

P°(cos 

♦  ) 

AXFN 

(N,  1, 

J)/AXFN 

(N,  2 

AXFN 

(N,  3, 

J) 

AXFN 

(N,  4, 

J) 
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The  functions  QZETA,  QCSPSI,  QSBF,  QALF,  which  generate,  respectively, 

£,  cos  \|(,  j^,(£)  and  P®(cos  ^ ) ,  are  also  called.  The  data  deck  for  XSCTN2 
contains  first  various  constants  AXFN  for  the  translational  integrals. 

These  are  followed  by  a  block  of  quantities  that  are  either  fixed  or  have 
proven  relatively  standard: 

NGLQ :  Number  of  Gauss -Laguerre  pivots  (GLQPVT)  and  weight  (GLQWT)  to 

be  subsequently  read  in. 

EMCRT:  Criteria  applied  to  error  measure  in  QGLQ. 

GLQPVT,  GLQWT:  Pivots  and  Weights  for  QGLQ. 

ZSBF :  Value  of  argument  modules  for  QSBF  below  which  Taylor  series 

is  used . 

SBFCV:  Criteria  applied  to  convergence  estimate  for  series  in  QSBF. 

CNCLT:  Criteria  applied  to  cancellation  estimate  for  explicit  evalu¬ 
ation  in  QSBF. 

E:  Criteria  applied  to  truncation  error  estimate  in  INTGRT. 

NGLM:  Number  of  modified  Gauss -Laguerre  p ivots  (BLMPVT)  and  weights 

(GLMEWT)  to  be  subsequently  read  in. 

NSTAT :  The  (odd)  number  of  reporting  stations  desired  on  the  trajectory. 

MODE:  (=1)  no  Rosen-Zener  model  computed,  (=2)  Rosen-Zener  model 

computed,  (=3)  as  2  but  intermediate  output  (associated  with 
NSTAT)  omitted. 

NEXT:  (=1)  a  news  ve loc ity /energy  (V,KEV)  read,  (=2)  a  new  "reaction 

block"read . 

NXCH:  (=1)  bypasses  proposed  "exchange"  modification. 

NELCN:  (=1)  elastic  matrix  element  correction  not  performed,  («2 )  cor¬ 
rection  performed. 

Following  these  quantities  are  one  or  more  reaction  blocks  and  one  or 
more  velocities  or  energies  depending  on  value  given  to  NEXT. 

APPN:  (-0)  causes  immediate  termination  of  execution,  (“X,  where  X  is 

any  alphameric  symbol  except  0)  causes  output  to  be  identified 
as  Xs  approximation. 

TAR:  Alphameric  identification  of  target  species. 

ETAR:  See  C  statement  following  read. 

DUM:  Dummy  variable. 

NTEL:  Number  of  active  target  electrons  (unused). 

ZTAR:  Net  charge  in  a.u.  on  target  after  collision. 

T0,VT:  Target  orbital  and  potential  parameters  (see  reference  1) 

PRO:  Alphameric  identification  of  projectile. 

EPRO:  See  C  statement  following  read. 

PRMSS :  See  C  statement  following  read. 

NPEL:  Number  of  active  projectile  electrons  (unused). 

ZPRO:  New  charge  in  a.u.  on  projectile  before  collision. 

P0,VT:  Projectile  orbital  and  potential  parameters. 

RHOA,  RHOB:  pa  and  pj,  referred  to  above.  If  RHGA=0,  RHOB  is  assumed 
to  contain  y. 

EFFRAD :  Assumed  maximum  distance  at  which  systems  interact. 

V.KEV:  Velocity  in  a.u.  and  energy  in  KeV  of  collision,  one  should  be 
zero. 

The  output  of  the  program  is  labeled  so  as  to  be  self-explanatory  on 
cards  and  printer  both. 
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FLOW  CHART  FOR  XSCTN2  MAIN  PROGRAM 


Read  quantities  associated  with  nontranslational  integrals. 

Read  quantities  associated  with  translational  integrals. 

Read  operational  parameters:  NSTAT,  MODE,  NEXT,  NXCH,  NELCN 
■►Read  first  part  of  "reaction  block"  data:  parameters  of  atomic  systems. 

Test  APPN  ZerV  exit. 

^Nonzero 

-►Read  second  part  of  "reaction  block"  data:  collision  paraneters . 

Compute  V  or  KEV  given  one. 

Write  reaction  specification. 


Test  RHOA - - 

Nonzero  | 

Call  MLTI04  twice  generating 
CA,  CB  for  RHOA,  RHOB.  Compute 
GMM.  | 


Zero 


Interpret  RHOB  as  GMM. 


Set  up  pivoted  impact  parameter  VRHO  appropriate  to  GMM. 

Call  MLTIOA  repeatedly,  computing  numerical  (CPR) ,  Rapp-Francis  (RAPP) 
and  model  (ANA)  capture  probabilities  for  approximately  first  5/8  of  VRHO. 
Write  final  output  for  reaction  consisting  of:  reaction  specification; 
RHOA,  RHOB,  CA,  CB  and  GMM;  a  table  of  impact  parameters,  3  capture 
probabilities  and  their  contributions  to  the  total  cross  sections;  final 
statement  of  3  approximations  to  the  cross  section. 

Teat  NEXT . 
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FLOW  CHART  FOR  ML.TI04 


Enter  MLTIOA 
Is  RHO  >  EFFRAD? 

I 

Write  specification  of  trajectory.  Write  specification  of  trajectory 

Compute  half  trajectory  length  and  reason  for  skipping,  exit. 

XF  from  RHO  and  EFFRAD.  Set  up  _  ^ 

Return. 

reporting  stations  STAT  along 
trajectory.  Initialize  solution 


YO,  coordinate  XO,  Interval  HO. 
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FLOW  CHART  FOR  MLTI04  (cont.) 


Compute  CFR 
Test  MODE 

^  U,3 

Construct  fitting  parameters  kj  (KAY1)  and  k2  (KAY2 ) 

Compute  corresponding  analytic  solution  and  asymptotic  value. 

* 

Compute  RAPP 
Test  MODE 


k: 


Write  basic  Information  about  calculation  along  trajectory  f  on  T 
Test  MODE 


Write  additional  information  regarding  modeling,  compute  ANA. 

I 

Write  modeling  parameters. 

Write  3  versions  of  calculated  capture  probability  CPR,  RAPP,  ANA 


Return. 


APPENDIX  II 


STRIPPING  CROSS  SECTION  PROGRAM 


The  program  that  calculates  the  electron  stripping  cross  section  con¬ 
sists  of  a  very  short  driver  and  six  subprograms.  The  program  Is  written 
for  the  CDC6000  series  computers  In  FORTRAN.  Actual  testing  and  production 
were  performed  on  a  CDC6400  machine.  The  calculations  required  little  com¬ 
puter  time,  only  about  30  seconds  being  used  to  calculate  the  stripping 
and  energy  loss  cross  section  for  about  10  different  projectile  systems  in 
one  target  gas  system  for  a  wide  range  of  impact  energies.  The  program  is 
simple  and  straightforward,  and  contains  numerous  comment  cards,  so  that 
exceptionally  detailed  program  documentation  is  not  required.  A  schematic 
flow  diagram  of  the  program  is  shown  in  Figure  67. 

The  main  program  STIP  calls  the  important  working  subroutines  and  reads 
the  parameter  N, which  is  on  the  final  data  card  for  each  calculation.  The 
value  to  which  N  is  set  controls  the  reading  of  the  next  data  set.  All 
new  data,  new  target  data  only,  or  new  projectile  data  only  Is  called  for, 
with  N  ”  1,  2,  or  3,respectively.  If  N  is  zero  or  negative  the  program 
stops.  Subroutine  REDD  reads  the  input  data  under  the  control  of  the 
parameter  N.  Subroutine  WORK  calculates  the  limits  of  integration  and 
computes  the  electron  stripping  and  energy  loss  cross  sections  by  calcu¬ 
lating  the  values  of  the  integrals  using  Gauss  gradratures.  Subroutine 
WORK  also  controls  the  printing  of  the  results. 

Subroutine  TERP  prepares  the  interpolation  arrays  EFF  and  SEFF  for  the 
elastic  and  Incoherent  scattering  factors  from  the  input  arrays.  Subroutine 
FIT  actually  calculates  the  array  elements  for  the  interpolation  arrays. 

A  cubic  fit  is  used  for  the  interpolation.  Function  S  calculates  the  in¬ 
coherent  scattering  factor  for  a  given  momentum  transfer  by  using  the  inter¬ 
polation  array  SEFF  for  the  cubic  interpolation.  Function  F  calculates 
the  elastic  form  factor  using  the  array  EFF. 

Hie  major  variables  in  the  various  segments  of  the  program  are  identi¬ 
fied  as  follows; 

N  is  the  parameter  that  controls  repeated  executions  for  different 
target  and  projectile  systems  as  described  previously. 

NINT  is  the  number  of  points  used  in  the  Gauss  Integrations. 

WT(I)  is  the  I-th  weight. 

PIV(I)  is  the  I-th  Gauss  pivot. 

NPE  is  the  number  of  input  data  points  for  the  elastic  form  factor. 

ITYP1  and  ITYP2  are  parameters  which  allow  the  input  data  to  be  in  the 
different  units  commonly  used  for  the  scattering  factor  data.  De¬ 
tails  are  on  the  comment  cards  in  the  subroutine  REDD. 

XKE  is  the  one -dimensional  array  containing  the  tabulated  values  of  K 
for  which  the  elastic  form  factor  is  given. 
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XFKE  is  Che  array  of  elastic  form  factor  values. 

NPI  is  the  number  of  input  data  points  for  the  incoherent  scattering 
factor. 

XXI  is  the  array  of  K  values  for  which  Incoherent  scattering  factor 
data  are  given. 

XFKI  is  the  array  of  incoherent  scattering  factor  data. 

IZ  is  the  number  of  electrons  in  the  target  system. 

TMASS  is  the  mass  of  the  target  system  in  atomic  mass  units. 

EXCT  is  the  target  excitation  energy  as  discussed  in  Section  III-2. 

NSP  is  the  number  of  shells  in  the  projectile  system. 

NES  is  the  array  containing  the  numbers  of  electrons  in  each  projectile 
shell. 

POT  is  the  arrary  of  ionization  potentials  for  the  various  projectile 

shells . 

PMASS  is  the  projectile  mass  in  atomic  mass  units. 

EMIN  is  the  lowest  laboratory  collision  energy  in  electron  volts. 

EMAX  is  the  highest  laboratory  collision  energy  in  electron  volts. 

DELTAE  is  the  energy  increment  for  the  energy  range. 

RM  is  the  reduced  mass  of  the  target-projectile  system. 

ECM  is  the  center  of  mass  system  energy  in  electron  volts. 

EAU  is  the  laboratory  energy  in  atomic  units. 

V  is  the  nuclei  relative  velocity  in  atomic  units. 

ULE  is  the  upper  limit  of  Integration  for  the  elastic  scattering  con¬ 
tribution. 

ULI  is  the  upper  limit  of  integration  for  the  inelastic  scattering 
contribution. 

WIE  is  the  lower  limit  of  integration  for  the  elastic  scattering  con¬ 
tribution. 

WLI  is  the  lower  limit  of  integration  for  the  inelastic  scattering 
contribution. 

CROSS  is  the  array  that  holds  the  contributions  to  the  stripping  cross 
section  for  each  shell. 

CROSI  is  the  array  that  holds  the  contributions  to  the  energy  loss 
cross  section  for  each  shell. 

TXS  is  the  total  stripping  cross  section. 

TSY  is  the  total  energy  loss  in  stripping  cross  section. 

EFF  is  the  array  of  parameters  used  in  the  interpolation  for  the  elastic 
form  factor. 

SEFF  is  the  array  of  parameters  used  in  the  interpolation  for  the  inco¬ 
herent  scatter  ng  factor. 

Following  the  listing  of  the  program  and  subprograms,  a  listing  of  data 
for  the  calculation  of  the  stripping  and  energy  loss  in  stripping  cross 
sections  for  hydrogen  projectiles  in  argon  gas  is  given.  A  detailed  des¬ 
cription  is  not  presented  because  the  conment  cards  in  the  subroutine  REDD 
give  sufficient  information. 
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PROGRAM  ST  IP  1  INPUT .OUTPUT < TAPE1  1  I NPUT  *  T  APE  2  =  OUTPUT) 

't-  ThTTS  PROGRAM"  CALCULATES  THE  5TR  rPElTJC  CROSS  SEC  T  I  ON  US  TNG  A  MODIFICATION 
C  OF  THE  DMITRIEV-NIKOLAEV  FORMUL A T I ON . JE T P  17,447(1963) 

- Fm - — — ■  — -  - . — 

3  CALL  REDDIN) 

~C  N-l  TO  READ  All  NEW  DATA 
C  N«2  TO  READ  NEW  TARGET  DATA 
~C"N«3  TO  READ  NEW  "PRO  JETTl  LE  '  DATA 
C  WE  always  read  a  new  energy  range 


CALL  WORK 
~  READ  rr,  17  N 

1  FORMAT! 14) 

-  mw7“2rr,3 

2  STOP 


— 

SUBROUT  INF  REDD ( N 1 

COMMON/lNT/  W  T ( TOO  >  ,  P'1  V  Cl  00 1  »  N  I  NT 

COMMON/ TAR /WORD ST ( 20)  .NPE.XXEt 100)  »  XF  XE ( 100)  .NPI  .XXl  ( 1 00 >  » XF XI  (  100 

-  - 

C 

C 

11.12 .TMASS.tXCT 

COMMON /PRO J/WOR DSP ( 20 1 . NSP . N E S ( 2 0 1 . POT ( 2C ) .PMASS 
COMMON/ENGy/EMIn.EMAX.DELTAT 

THIS  SUBROUTINE  READS  THE  INPUT  DATA 

GO  TO  (1.2.31  N 

DATA  FOR  NUMERICAL  INTEGRATION 

-  -  - 

— 

c 

1  R  E  AD (1.41  NlNl 

A  FORMAT (  I  A ) 

NINT  IS  THE  NUMBER  OF  WETGHfS  AND  PIVOTS'  USED 

DO  5  1=1. NINT, 1 

R E AD  ( T, 6  1  WTIT)  .'Pi  V(  D 

5  CONTINUE 

- 

— 

c 

c 

c 

s  Format < 2 t E  1S.H , lx i 1 

WT  AND  PIV  ARE  THE  WEIGHTS  AND  PIVOTS  RESPECTIVELY 

READ  THE  TARGET  DATA 

2  RE  AD (1,71  WORDST 

7  FORMAT( 1 UA8 1 

WORDST  ARE  TWO  ALPHANUMERIC  CARDS  GIVING  DETAILS  OF  THE 

TARGET  SYSTEM 

c 

c 

c 

READ (1,8)  NPE , I TYP1 , I t YP2 

8  FORMAT ( 3 (  I  A, IX 1  1 

NPE  IS  THE  NUMBER  OF  DATA  POINTS  FOR  THE  ELASTIC  FORM 
ITYP1=1  FOR  X  IN  ATOMIC  UNITS. NUT  1  FOR  <  IN  TERMS  OF 
I TyP2=1  FOR  F(0)=1  ,  NOT  1  FOR  F(0)«Z 

READ ( 1 , A  1  IZ 

FACTOR 

S I N ( THETA) /LAMBDA 

- c 

12  IS  THE  NUMBER  OF  ELECTRONS  IN  THE  lARGET 

DO  9  1=1. NPE. 1 

READ!  1,6)  XXE ( 1  1  .XFxE (  I) 

9  CONTINUE 

I  F  (  ITYP1-D10.11.10 

10  DO  12  I-l.NPE.l 

1  2 

XXE ( I 1 -6.6A971*XXE ( I ) 

12  CONTINUE 

11  IF (  ITYP2— 1  1  13.1A, 13 

IT  DO  15  1=1 .NPE , 1 

XFXE( I >=XFKE( I 1 /FLOAT ( I 2 1 

15  CONTINUE 

8 

c 

c 

c 

1 A  CONTINUE 

R E AD (1,8)  NPI  .  I  TYP1 , I TYP2 

NPI  IS  THE  NUMBER  OF  DATA  POINTS  FOR  THE  INELASTIC  FORM  FACTOR 

ITYP1  =  1  FOR  K  IN  ATOMIC  UNITS. NOT  1  FOR  X  IN  TERMS  OF  S I N (  T HE T A  I / L«MbD A 
ITYP2  =  1  FOR  S I N C ( I  N  F  )  =  1  ,  NOT  1  FOR  SlNC(lNF)=2 

DO  16  1=1, NPI. 1 

5 

READ (1,6)  XX I  (  I  )  ,  XFX I  (  I  ) 

16  CONTINUE 

I  F (  ITYP1-1  H7.18.17 

17  DO  19  1=1 »NP I . 1 

XXI  <  I  )  * 6  •  6  A"97 1  * XX I  (  i  ) 

19  CONTINUE 

3 

18  IF ( I TYP2-1 ) 20,2 1 ,20 

20  DO  22  1=1 .NP I . 1 

xfx  i  ( i » =xfx  i  ( n /float < i2i 
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22  CONTINUE 
2 r READ (1.6)  TMASS 
C  TMASS  IS  THE  TARGET  MASS  IN  AMU 
L  U»T»  KtADlNu  SlAlLMtNTb  F  UK  I  AKUt  I  CUMMLt  I  tl> 

READ (1.6)  EXCT 

C  FXCT  IS  A  LOWE5T  EXCITATION  ENERGY  FOR  "THE- TARGET  SYSTEM 
3  GO  TO  ( 23 . 2 A . 23  )  N 
25  RE~AO  (T.7 )  W0RD5P 
C  WOROSP  ARE  TWO  ALPHANUMERIC  CARDS  GIVING  DETAILS  OF  THE  PROJECTILE  SYSTEM 

- READ  1  r»y>  NSP - 

C  NSP  IS  THE  NUMBER  OF  SHELLS  CONSIDERED  IN  THE  PROJECTILE 

DO  25  I-l.NSPn  -  -  -  —  “  -  - 

READ ( 1 .26 )  NES<  I  >  .POT <  1  ) 

25  CONTINUE  '  - 

26  FORMAT ( IA. 1X.E15.8 ) 

"CUPS'  T5'  THE  Number  OF  EL"ECTR0N5  IN^rWE-5HETn: - 

C  POT  IS  THE  THRESHOLD  ENERGY  TO  REMOVE  THE  ELECTRON. IN  AU 
READ (T. 6 )  PMASS 

C  PMASS  IS  THE  PROJECTILE  MASS.  IN  AMU 

_  2a  continue  —  ■ 

READ (1.27)  EMIN .EMAX .DELTAS 


27  F0RMAT(3(E15.8.iX) ) 

C  EMIN  IS  THE  LOWEST  COLLISION  ENERGY  CONSIDERED 
C  EmAx  IS  THE"  HIGHEST  COLL  I  ST  ON“ENERGY  COFTSIDETTETD 
C  DELTAS  IS  THE  INCREMENT  IN  COLLISION  ENERGY 
RETURN 
END 


SUBROUTINE  WORK 

COMMON/lNT/  WT  (  lOO )  ,PI V!  ICO) *  N I NT 

COMMON /PROJ/WOR DSP  ! 20  > . NSP . NE S I  20  >  .POT ( 20  )  .PMASS 
COMMON /EnGY /EMIN. EMAx. delta e 

COMMON/ TAR /WORDS T (  20)  .NPE.XKE(IOO)  , XFKE  (  1 oO )  .NPI  .XKI  (100 >  .XFKI  (lOO 
11.12 . TMASS. ExCT 
DIMENSION  CROSS ( 20  )  .CROSl  ( 20) 

C  THIS  ROUTINE  ^ALTuLAtES  THE  CROSS  srcnoN 
WRITEI2.il 

1  FORMAT! 1H1) 

WRITEI2.2) 

2  FORMAT! 22H  TARGET  f  DEN  T  IT  rCAfl  ONT 
WRITEI2.3)  WORDST 

3  Format (Tx . ioas ) 

WR ITE! 2 .A ) 

- - a  FORMAT  (  75R  "PT5CJECT  T  lE  IDENTIFICATION) - - — - - - - - - 

WR I T  E  <  2 ♦ 3 )  WORD  SP 
RM=PM ASS* TMASS/ ( PMASS+TMASSI 
C  EMI N .EMAX .DEL TAE  ARE  IN  EV  IN  THE  LABORATORY  SYSTEM 
E  =  EM  I N 

5  ECM=E*TMASS/(PMASS+TMASS ) 

EAU=E/?T.21 

C  ECM  IS  THE  CM  ENERGY  IN  EV.  EAU  IS  THE  LAB  ENERGY  IN  AU 
V=SQRT( 2.0*EAu/ (PM ASS* 18 22. 3  1  1 
C  V  IS  THE  RELATIVE  VELOCITY  .IN  Au 
ULE*2.0*V 
UL  I  *=  V 

- wirrmTTST  etfthvv - - - - 

6  FORMAT ( 1 1H  E(LAB,EV)  =  .E1S.8,10H  E  (  CM , E V )  =  . E 1 5  -  8 . 2 3H  RELATIVE  VELOC 
1 ITY! Au)=,El5.8 ) 

DO  7  1=1, NSP. 1 
WLE=S0RT  f2 • 0*P0  T < I ) ) 

IF  <  1.0-2.0*EXCT /(V*V)  )  100.101.101 

100  WLI'2.0*ULI  - 

GO  TO  10 

101  AA*V* ( 1 .O-SQRT!  1.0-2.0*EXCT/ ( V*V )  )  ) 

IF(AA-WLE)  8,8.9 

8  WLI=WLE 
GO  TO  10 

9  WLI*AA 

10  CONTINUE 

C  PERFORM  THE  ELASTIC  INTEGRATION 
IF  fULF-WLF 111,11.12 

1 1  CROSS (  I  1=0.0 
CROSl (  I  ) =  0 . 0 
GO  TO  7 

12  SUM 1=0.0 
SUM3=0.0 

DO  13  J=l.NlNT.l 

X  =  0 • 5* ( ULE  +WLE ) +0.5*rULE-WLEl*PlV( Jl 
22 Z  =  F( X  ) 

2  2  =WT ( J ) * (  (  1.0-ZZ?  ) **  2 )  /  (X*X*X) 

221=0. 5*ZZ*X*X 
SUM1 =SUM1+ZZ 
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_  SUM3  *SUM3+ZZ ^  _ 

13  Continue 

SUM1«SUM1*0.5*!ULE-WLEI*FLOATI i  z*  i  z i 

- 5UHygp'^TStjH3jHULE.jWLb)«FLOAT(  1Z»1Z) - 

C  DO  THE  INELASTIC  INTEGRATION 

“TnUtT-WLir  T4.T4.T5 - 

14  SUM2-0.0 

- 5U«4^^.0 

GO  TO  16 

- 15  5UM2-0.0 - 

SUM4-0.0 

oo  i7  jir.wrwr.T~  - - 

X«0.5*iULl+WLl  )+0.5«<ULl-WLl  9*PIvl J> 

rzi^si xf 

ZZ-WT! J)*ZZZ/IX*X*X) 

- rzzi-a.~5»zz'rx*r - 

SUM2-SUM2+ZZ 

- 5UH4“ SUM  4 ♦ ZZ  2 - 

17  CONTINUE 

- 5UM2-5UMZ»  0V5F(  ULT-~WLT7«F  LUA  T  ?  ITT - 

SUM4«0.5*SUM4*< ULl -WLI  I*FLOAT(  IZ  > 

16  SUM-SUH1+SUH2 

SUMP-SUM3+SUM4 

_  CROST  (  r7-5UMP*7  .~rr3773E-I6»FLOATrNEr5TT7  1/!V*V) 

CROSS!  I  )»SUM*7.03773E-16*FLOATINES( I  I  I / ! V*V 9 

- -  •  T CONTINUE - 

TXS-0.0 

- T5Y5O20 - - — - - 

DO  18  I-l.NSP.l 

- - T5Y*T  5Y-FCR05T  ("IT  - 

TXS»TXS-*-CROSS! I ) 

- IT  CONTINUE - 

WR  I  T E ( 2 . 1 9 )  TXS 

- T9  FOHHA'T  1 24H  IOTAL  IROSVSECTIUN  TS  .'E13.B1 - 

WR I T E I  2 • 24  )  TSV 

-  74  FORMAT  (7aH~~ENEROVCO~5S~  CRUS5  ~SECT  TOW  ~TS  .El  5. 8) 

WR  ITE ( 2 .21 ) 

- - fnrjIXrTiZUTT  CROS5ITr7I»TVNSP9  - 

WR I TE ( 2 • 21 ) 

- WRTTFT7T7U7  rCNOSTTTTTTWTTlfSPl - 

20  FORMAT (6! 1X.E15.8)  ) 

21  FORMAT !28H  CONTRrBUTTONS- OT  THE  SHELVS7 
WR ITE! 2 .22  9 

~27  FORMAT  TIFT  7 - - 

E-E+OELTAE 

- rriE-EMAX  9  ~S,5  .'73 - — - 

23  RETURN 
~  ~  ~  EnD 
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SUBROUTINE  TERR _ 

COMMON  /  T  AR /WOROS  T ( 2 0  I  . NPE • XKE (  100)  .XFK£"<  100)  .NP1  .XKI  (  100)  ,XFKI  (  100 
1  )  < IZ.TMASS.EXCT 

COMMON/ TP T /EFF (4,100).SEFF(4,100) 

COMMON /EZP  /CEX_TP 

C  THIS  ROUTINE  CALCULATES  THE  INTERPOLATION  'ARRAYS 
C  FOR  THE  ELASTIC  FORM  FACTOR 
MAX«NPE-3 
DO  1  1=1. MAX. 1 

- X  1 1  =XXE  <  I  1 - 

X  2 1 =X  KE (  1  +  1  ) 

X  3  1  =  X  K  E  (  1  *  2  > 

X41*XKE^I  +  3  )_ 

F 1-XFXE (  FI 
F2-XFKE (  I  +  1  ) 

- F3=xP'kE< 1+5  i - 

F4=XFKE< 1+3 ) 

■  CALL"  F  I  T  (  X  1 1  .  X 2  1  ,  X  3 1  ,  X  A  1  .  F  r, F 2  ,T 3  . F  4  ,T5 1  ,  D 2  ,  D 3  »  D4  ) 

EFFll.I ) =D 1 
EFF  (  2 . D  =D2 
EFF ( 3 . I )=D3 
EFF) 4, 1 ) =04 

1  CONTINUE 

c  calculate  the  Extrapolation  constant  cextp 

CEXTP=(-XFK.E(NPE  )+SORT  {  XFKE  (  NPE  )  )  )  /  (  XFKE  (  NPE  )  *XKE  1  NPE  )  *X  XE  (  NPE  )  ) 

c  For  the  Inelastic  form  factor 

MAX  *NP 1  —  3 
DO  2  1*1 .MAX . 1 
X11=XXI ( I »_ 

X  2  1  *  X  X  I  (  1  +  1 ) 

X31»XXI  I  1+2  ) 

X4l=XKI  (  1  +  3  ) 

Fl-XFKI ( I ) 

- F__^_TTT_n - 

F3=XFXI ( 1+2 ) 

F4*XFK I ( 1+3 ) 

CALL  FIT(X11.X21.X31.X41.F1.F2.F3,FA,D1.D2,D3.D4) 

SEFF( 1 . I ) *D1 
SEFF (2,  I  ) =02 

SEFF (3,1) =  D  3  "  ~  -  '  ’  - 

SEFF (A. I >=D4 

2  CONTINUE 
RETURN 
END 


SUBROUTINE  F I T I  X 1 1 « X 2 1 • X 3 1 . X4 1 . F 1 . F 2 , F 3  ,  F4 , D !  •D2,D3,D4> 
X 12-X1 1*X  1 1  " 

X 1 3*X  1 1*X 1 2 


?*XZ1*X21 
X23*X21*X22 
X32«Tt31*>T3  1 
X33-X32*X3  1 
X42*X4l*X4  I 
X43-X41*X42 

02-<X12-X22'mxrr-X'yTT-m2~X32  )*(  Xll  1 
03>^(XI3-X23  >*<X  11-X31  >-(X13-X33)*(XlT*‘- 

04  •  fp y-f ?  )  *  ( x  n-x3i )  -  it  r-rfiFiyri-K? 1 

S2*  <  X 12-X  32 >*<Xll-X4l»-<X12-X42)*(X:t-»  ,i 
53*  fX  T3^  >T3  3  I »  EX  Tl  -  X4  T >  - 1 X  TT-  X4T) »  rx  '  l  j  ;  ) 

S4 * l F 1— F 31*1X11— X41  1-IF1-F4  >*<X11-X3m 

154-  (04*  52-54*02  J 77 03»,52-57*02l -  — 

D3*(Q4-D4 *031/02 

02*  TF1-F  2-04*  (XT3-X25I  -03»fX  T2-X2 T\ T/1  SO  1  -X  21 1 
D1 *F 1— D4*X 1 3— D3*X 1 2  — 02*X 1 1 

RETURN  - - - 

END 


FUNCTION  SIX) 

COMMON/ T  AR /WORD ST (20) .NPE.XkE(  100)  .XKk£ (  i 00  )  .NPI  ,XKI  I  100 )  ,XFKI  (  100 
1 ) , 12 .TMASS.tXCT 

COMMON  /  Tt>T  /EFF  <4,100).SEFF(4.100) 

COMMON/EZP /CEXTP 

C  THIS  FUNCTION  CALCULATES  THE  TNELASTTC  FORM  FACTOR  FROM  THE  ARRAY-  5EFF 
I  F ( X )  1.2*3 

1  WRITE!  2.4") 

STOP 

- 4  FORmA  t  T53H  NFGXTTvE'  MOWFflTOH  TRANSFER  ERROR! - 

2  IF (XK I ( 1 ) )  5,5.6 

5  S*XFK  I  <1 ) 

RETURN 

6  S  =  SEFF'(T.1>  +  SEFF  I  2  ,  n*X  +  SCFF(  T.Tr*)C**X+SEFF  (  4, 1  )*X*X»X 
RETURN 

3  I F  (X-XK  1(1)1  6 . 5  _____  _____  -  —  —  - 

7  I  F (X-XK 1(2))  6 ,8 ,9 
8“  S  =  XFX  1(2) 

RETURN 

9  I  F (X-XKl (  3  )  )  6.  TO, 1 1 

10  S=XFkI ( 3 ) 

EFTUEn 

11  IFIX-XKl (NPl ) )  12.13,14 

14  FE2P *1.0/1 ( 1.0+CEXTP*X*X)**2 ) 

S* 1 «  0- < FE2P*F  E  2  P I 

RETURN 

13  S«XFKI  (NPI  ) 

p  p  turn 

12  MAX  *NP 1—3 

IF  (X-XK  I  (NPr-D  )  15,16.17 

17  S=SEFF ( 1 .MAX ) +  SEFF ( 2, MAX )*X+SEFF (3,MAX)*X*X+SEFF (4 .MAX ) *X*X*X 
RETURN 

16  S*XFK I (NPI -1 > 

—  RETURN 

15  IFIX-XKI (NPI— 2) )  18,19.17 

19  S=XFKI (NPI -2) 

RETURN 

18  1=2 

I  I  =4 

- 26  I F  ( X -"xK Till  IT ~TU  .21  .22  ' — - ' — " ' ' -  - —  - 

20  S*SEFF  (1,1)  +SEFF (  2,1)  +X+SEFF  (3.I)*X*X+SEFF(4,I)*X*X*X 
RETURN 

21  S=XFKI (II) 

RETURN 

22  I F (  I —MAX  123.24,24 

24  WRITE! 2 ,25 ) 

STOP 

25  FORMAT! 17H  LOGIC  ERROR  T N  S) 

23  1=1+1 

I  1  =  1  1  +  1 
GO  TO  26 


FUNCTION  F(X> 

COMMON/ TPT/EFEi  A.  ICO  I  .5EFF  (  4 , 100'T 

COMMON /TAR /WORD ST ( 20  > ,NPE ,XKE <  1  00  >  .  XF ICE  <  1  00  >  . NP  1  . X X  I  (  1  00  >  .  XF K  I  (  1 00 

- I~T  .  T7~<Tffir5T»EyCT - — 

COMMON / EZP /CEX  TP 

C  THIS  FUNCTION  CALCULATES  THE  EL  AST  ITT  FORM-  FACTOR  FROM  THE  ARRAY  EFT 
IF (X  1  1,2,3 
1  WRITE (2, Al 
STOP 


*  f  'li  ATJ2M|  *  r»]  •  i  ai  illuHIC*#' Viw139IiBSI'lk29]*SJ 


2  IF (XKF (  1  )  I  5,5,6 

5  F-XFKE(l)  ’  -  - 

RETURN 

6  F-EFFTT.  II  -t-EFFI  7,1  I  *X+EFFT3i  1  )»X»XTEFF  T4, 1 ) *X*X*X 
RETURN 

~5  Tmr=yiErm i  s.r.7 - 

7  IF (X-XXE (2)1  6 .8 .9 
ft ‘  F«XFXE<2  ) 

RETURN 

9  TFIX-XKEI3I )  6.10,11 

10  F«XFKE<3) 

RETURN  ‘ 

11  IF(X-XKE(NPEI )  12,13,14 

14  F  * 1.0/1  (  1.0*CEXTP*X*X)»»2) 

RETURN 

13  F4TCFKF(NPE) 

RETURN 

12  MAX  «NPt  — 3  - - — _______  - - -  -  - 

I  F ( X-XKE (NPE-1 1  )  15,16,17 

17  F-EFFT  I  ,MAXI*EFF(2,MAX  J«XTETF(3,PAX  J*X*X  +  EFF  (4,  MAX  )  *X*X*X 
RETURN 

16  F«XFKE(HPE-ri 
RETURN 

15  IF  IX-XKE  INPE-21  )  18.19.17  - - — - - - - - - — 

19  F«XFKE (NPE-2) 

RETURN 

18  I  «  2 

I  I  -4 

26  IF (X-XKE < 1  I > )  20,21,22 

20  F«EFF  (  1  ,  I  IFEPFI  2,1  )  # X  +E FF' (  3  .  IT  *X * X  +F FF  (  4  ,  TT»5T*'JT»x—  — 
RETURN 

21  F«XFXE(  I  I  ) 

RETURN 

22  IF( I-MAX )23,24.?4 

24  WR  ITE ( 2 .25  I 

- STOP - 

25  FORMAT!  17H  LOGIC  ERROR  IN  FI 

23  I  -  I ♦ 1 

I  l  «  I  I  +  1 
GO  TO  26 
END 
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